Spectral reduction: a statistical description of turbulence 
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Abstract 



A method is described for predicting statistical properties of turbulence. 
Collections of Fourier amplitudes are represented by nonuniformly spaced 
modes with enhanced coupling coefficients. The statistics of the full dynam- 
ics can be recovered from the time-averaged predictions of the reduced model. 
A Liouville theorem leads to inviscid equipartition solutions. Excellent agree- 
ment is obtained with two-dimensional forced-dissipative pseudospectral sim- 
ulations. For the two-dimensional enstrophy cascade, logarithmic corrections 
to the high-order structure functions are observed. 
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Many practical applications for spectral simulations of turbulence exist where it would 
be desirable to evolve modes that are distributed nonuniformly in Fourier space, devoting 
most of the computational resources to the length scales of greatest physical interest. This 
idea has led to the development of a new reduced statistical description of turbulence, called 



retained in simulations of turbulent phenomena. It exploits the fact that statistical moments 
are much smoother functions of wave number than are the underlying stochastic amplitudes. 

The concept of wave-number reduction is not new. In the method of constrained decima- 
tion [@-§]], a stochastic forcing term is added to model the effect of the deleted modes on the 
retained modes. She and Jackson have proposed a reduction scheme in which the linear (vis- 
cous) term is modified ||. In spectral reduction, a third alternative is chosen: the nonlinear 
coefficients are enhanced to account for the effect of the discarded modes on the explicitly 
evolved modes. There have been other more heuristic attempts at wave-number reduction 
0-0; these methods typically neglect nonlocal wave- number triad interactions (which play 
a particularly important role in two-dimensional turbulence). Unlike the renormalization 
group [TIJ method, which retains only large-scale modes and attempts to express the effect 
of the small-scale modes using a self-similarity Ansatz, spectral reduction retains certain 
modes from all scales, while discarding other modes associated with these same scales. The 
generality of the formulation allows one to refine the partition wherever the physics dictates. 

In this Letter we restrict our attention to homogeneous and isotropic incompressible 
turbulence in two dimensions. The appropriate spectral transform in this limit is the inte- 
gral Fourier transform, under which the two-dimensional Navier-Stokes vorticity equation 
becomes 



Here vu models time-independent linear dissipation or forcing and the interaction coefficient 
efcpq = (z-pXq) S(k + p + q) is antisymmetric under permutation of any two indices (* 
denotes complex conjugation and = indicates a definition). We restrict the integration to 



spectral reduction [[TJ, which dramatically reduces the number of spectral modes that must be 
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a bounded wave-number domain T> that excludes a neighborhood of k = 0. As a conse- 
quence of the antisymmetry of €k pg , in the inviscid limit [v^ = 0) Eq. ([]]) conserves the 
energy | f v dk \iOk\ 2 /k 2 and enstrophy | J^rffc \uik\ 2 . It is believed that the energy and en- 
strophy play fundamental roles in the dynamics of the turbulent cascade. 

We introduce an arbitrary coarse-grained grid that partitions T> into connected regions 
called bins. The bins are labeled by capital letters to distinguish them from the continuum 
wave numbers, which we represent by lower-case letters. To this grid, we associate new 
variables VLk = A^ 1 J A uj^ dk, where is the area of bin K. The exact evolution of VLk 
is given by 

^ + (^)k = E a ^(^k) > ( 2 ) 

Ul P,Q \ q I KPQ 

where {■) K denotes a bin average and the operator 

if) kpq = A A A / dk [ d P [ d <lf> ( 3 ) 
AkA p Aq J Ak J Ap J Aq 

depends only on the bin geometry. The geometric factors (f) kpq can De efficiently computed 
using a combination of analytical and numerical methods [p I |L3H . Since they are independent 
of both time and initial conditions, they need only be computed once for each new wave- 
number partition. The reality condition Qk = ^1k> where —K denotes the inversion of 
bin K through the origin, will be respected for partitions that possess inversion symmetry. 

Equation (fj) is unfortunately not closed. If Uk were naively approximated by its bin- 
averaged value Qk, one would obtain 

^ + <^^ = £ApA Q /^\ n* P n* Q . (4) 
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In the inviscid limit, Eq. conserves the coarse-grained enstrophy \^2, k \^k\ 1 ^k 
since (ek pq /q 2 ) kpq * s antisymmetric in K <-> P. However, the coarse-grained en- 
ergy \ J2k \^k\ 2 Ak/K 2 is not conserved since (^k Pq / q 2 ) KP q / K 2 is not antisymmetric 
in K <-> Q (here K denotes the magnitude of some characteristic wave number in bin 
K). However, both of these desired symmetries can be reinstated by replacing the factor 



i e kpq/ q 2 ) KP q in Eq. (|J) with the slightly modified coefficient {^k Pq ) KP q /Q 2 ■ The relative 
error introduced by this modification is negligible in the limit of small bin size, being on the 
order of the squared relative variation in the wavenumber magnitude over a bin. The result, 

^ + W k »k = £ ApA q ^^«^, (5) 
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which we call the spectrally reduced Navier-Stokes equation, is a more acceptable alternative 
than Eq. @ as a closure of Eq. @: not only does it reduce to the Navier-Stokes equation 
in this limit, but it also conserves both energy and enstrophy, even when the bins are large. 
The final modification leading to Eq. (|^) partially compensates for the error introduced by 
the crude approximation uj k ~ Qk and leads to the same general structure and symmetries 
as Eq. ([[]); in this sense spectral reduction may be regarded as a renormalization. 

If the bins are large, the true vorticity will vary rapidly with wave number within each 
bin and it is unlikely that Eq. ([|) will yield a reasonable description of the instantaneous 
dynamics. However, the time-averaged (or ensemble-averaged) moments of Eq. ([5]) satisfy 
equations that closely approximate the equations governing the exact bin-averaged statistics. 
For example, a time average (denoted by an over-bar) of the bin-averaged enstrophy equation 
derived from Eq. (|l|) leads to 
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If the true vorticity is a continuous function of wave number, the mean value theorem 
for integrals guarantees the existence of a wave number k in bin K such that Qk = w K . 



Furthermore, time-averaged quantities such as \ujk\ are generally smooth functions of the 
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wave number k. We thus deduce that \Qk\ — \^k\ ~ l^fel f° r an & in bin K. Similarly, 



the triplet correlation u^UpU* is a smooth function of k, p, q when restricted to the surface 
defined by the triad condition k + p + q = 0. 

To good accuracy the statistical averages in Eq. (^) may therefore be evaluated at the 
characteristic wave numbers K, P, Q of each bin, yielding 
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Moreover, to the extent that the wave-number magnitudes vary slowly over a bin, Eq. 
may equally well be reduced to the (nonlinearly conservative) approximation 
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which is precisely the evolution equation for the time-averaged enstrophy obtained 
from Eq. fl5|). Similar arguments for the higher-order statistical moments can also be made, 
suggesting that spectral reduction can indeed provide an accurate statistical description of 
turbulence, even when each bin contains many statistically independent modes. As the 
partition is refined, one expects the solutions of Eq. (|j) to converge to the those of Eq. (^j). 

In the absence of forcing and dissipation, the (untruncated) two-dimensional Euler equa- 
tions can be written in a noncanonical Hamiltonian framework jTJ| as u k = J dq J kq 5H/5u q , 



where H = ~ J dk \u k \ 2 /k 2 and J kq = j dpe kpq uj*. The Liouville theorem 
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then follows immediately from the properties of e kpq . 

When v k = 0, Eq. (|) can be written in a similar form as Qk = J2q Jkq9H / 'OVLq, 
where H = \Y,k \^k\ 2 ^k/K 2 and 

Jkq = ^ A P (e kpq ) KPQ fip. (9) 
p 

It is an open question whether (an untruncated version of) Eq. (£5|) satisfies the Jacobi 
identity, which would make spectral reduction a Hamiltonian approximation. What is certain 
is that the respective Liouville theorem 

^8QK = ^fdJKQdH_ d 2 H \_ 

K dnK hi ^ dnK dnQ KQdn Kdn Q ) 



is obeyed, as a consequence of the antisymmetry of (^k P q) KP Q- If the dynamics is mixing, 



the inviscid system will then evolve toward equipartition [|T|,|T6|]; this was verified for spectral 
reduction numerically, using a fifth-order conservative Runga-Kutta integration algorithm 
that conserves quadratic invariants to all orders in the time step [f[7],|]l|. When using nonuni- 
form bins it is necessary to rescale the time derivative d/dt in Eq. ([5]) to (A /Aj^)d/dt, where 
Ao is the minimum bin area, to obtain an equipartition of modal (instead of bin-averaged) 
energies. This modification to the transient evolution will be discussed further in a future 
paper. 

Upon adding to Eq. (H) a random stirring force for k 6 [5, 7] and adopting the dissipation 
function u k = v L 6{7 — k) + v u k 2 (Q is the Heaviside function), we graph in Fig. || the 
time-averaged saturated energy spectra for four wave-number partitions to test how rapidly 
spectral reduction converges. The excellent agreement demonstrated between the predictions 
of spectral reduction and a (computationally more expensive) full dealiased pseudospectral 
simulation is obtained without fitting or the introduction of adjustable parameters. We 
also show the predictions of the realizable test field model (RTFM) |TB| , using wavenumber 
binning and setting the eddy-damping multiplier in this heuristic statistical closure to one. 

High-order moments are also accurately described by spectral reduction. A quantity 
of interest is the angular average S n (r) of the n-th (time-averaged) moment of velocity 



increments \v(r) — v(0)\ n , or structure function. In Fig. ^| we illustrate the scaling with 
distance r of a typical high-order structure function, S w (r), for the runs depicted in Fig. [I]. 
Slight variations in the predicted large-scale velocities are evident as overall vertical offsets. 
Note that the (discrete) pseudospectral calculation is an approximation to Eq. ([[]) at the 
large scales. 

One can readily investigate high- Reynolds number turbulence with spectral reduction, us- 
ing a polar partition in which the bins are logarithmically spaced in the radial wave number. 
A saturated turbulent state can be evolved for thousands of eddy turnover times to obtain 
statistically meaningful moments for comparison with theoretical predictions. For example, 
Kolmogorov's idea of self-similar energy transfer in the inertial range |19| led Kraichnan 
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20| to propose a logarithmically corrected asymptotic form for the energy spectrum E(k) of 
the enstrophy cascade. In a simulation with viscous dissipation active only at the smallest 
scales (to yield a pristine inertial range) and forcing via a linear instability (negative Vk), 
we apply spectral reduction to demonstrate the recent extension E{k) ~ k~ 3 x~ 1 ^ 3 (k) of 
Kraichnan's result to the entire inertial range, where ki is the smallest inertial-range wave 
number, %(&;) = \ia(k/ki)+xi, and the positive constant xi is set by the large-scale dynamics 



2i| . We verify in Fig. ||| the linear behavior of [k E{k)\ with respect to \n(k/ki) } using 
the values k\ = 16.4 and xi = 0-67 determined by a least-squares fit. The inertial-range 
energy spectrum is thus well described by Kraichnan's logarithmically corrected k~ 3 law. 

For the second simulation, we demonstrate in Fig. |4] the linear behavior of [r~ n S n (r)] 3 ^ n 
with ln(ri/r) on the interval 0.043 < r < r\ = 0.26 for various values of n. The implied 
scaling S n (r) ~ r™[ln(ri/r) + Xn]"^ 3 > where x' n ls a constant, is in agreement with both 



the asymptotic theory of Falkovich and Lebedev p2| and the recent experimental results of 
Paret et al. ||23[| , lending support to the claim that there are no high-order intermittency 
corrections in two-dimensional turbulence. The universality of this result will be investigated 
in a future paper. 

In this Letter, we propose a new technique that dramatically decreases the number of 
degrees of freedom required to simulate homogeneous turbulence. The statistically stationary 
state described by Fig. [| which would require 2048 x 2048 dealiased (3071 x 3071 total) 
pseudospectral modes, can be successfully modeled using only 32 x 8 bins. A notable feature 
of spectral reduction that distinguishes it from other statistical theories of turbulence is the 
existence of a control parameter (bin size) that can be varied to increase the accuracy of a 
solution. Moreover, spectral reduction does not make a closure assumption on the triplet 



correlation Q* k Q* p Qq appearing in Eq. (|8]); it circumvents the closure problem entirely by 
reducing the number of triplet correlations to a tractable number, instead of eliminating them 
in favor of lower-order statistical variables. Unlike statistical closures, spectral reduction thus 
does not destroy the phase information embodied in the triplet correlation. 

Spectral reduction appears to be a promising candidate as a statistical description of 



turbulence. We propose that it could be used to assess the effect of various dissipation 
mechanisms in large-eddy simulations, as a subgrid model, or even as a substitute for full 
simulation of high-Reynolds number turbulence. However, as it does not provide explicit 
insight into underlying dynamical processes, spectral reduction should be considered more 
as a computational tool than as a true analytical theory of turbulence. The latter challenge 
still awaits us. 

We would like to thank John A. Krommes, for the idea of anisotropic bin-averaging 
of statistical closure equations, R. E. Waltz, for the suggestion of applying that technique 
directly to the Navier-Stokes equations, and C. E. Leith for discussions on the Liouville 
theorem. This work was supported by U.S. Department of Energy Contract No. DE-FG03- 
96ER-54346 and the Natural Sciences and Engineering Research Council of Canada. 
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FIGURES 
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FIG. 1. Comparison of the turbulent energy spectra obtained with 16 x 8, 32 x 8, 64 x 8, and 
16 x 16 (logarithmically spaced radial x uniformly spaced angular) wave-number partitions, the 
RTFM, and a full 683 x 683 dealiased pseudospectral simulation (1024 x 1024 total modes). 
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FIG. 2. Angle-averaged structure function Sio(r). 
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FIG. 3. Linearity of [k 3 E(k)]~ 3 with respect to ln(k/ki) for an enstrophy inertial range between 
k\ = 16.4 and k = 330. The solid triangles are the predictions of spectral reduction. 
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FIG. 4. Linearity of [r~ n S{r)f/ n with respect to ln(n/r) for 0.043 < r < n = 0.26. 
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